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Forewarci 

This report, prepared by the Structural Dynamics Group, 
Department of Civil Engineer iing, College of Engineering and 
Applied Science, University of Colorado under contract NAS8-31619, 
presents the results of a study that investigated the conceptual 
approach termed ’’Coupled Base Motion Response Analysis of 
Payload Structural Systems". The study was performed from 
July 1975 through March 1976 and is presented in three sections. 

Section 1. Analytical Developments 

Section 2. Example problem typifying HEAO-MECO event 

Section 3. Program Code 

The study was administered by the National Aeronautics and 
Space Administration's George C, Marshall Space Flight Center, 
Huntsville, Alabama, under the technical direction of Mr. Wayne 
Holland. 
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Principal Nomenclature 

Support system 

Branch, system 

Primed branch system 

Masses of the support i^ystem 

Masses of the branch system 

Mass of the cantilevered branch system 

Stiffness of support spring 

Stiffness of branch spring 

Stiffness of primed branch spring 

Viscous damping coefficient for support system 

Viscous damping coefficient for branch system 

Viscous damping coefficient for primed branch system 

Damping ratio for support system 

Damping ratio for branch system 

Damping ratio for primed branch system 

Circular frequency of support system 

Circular frequency of branch system 

Circular frequency of primed branch system 

Reduced mode shapes of support system with respect 
to Interface 

Mode shapes of cantilevered branch system 
Mode shapes of primed cantilevered branch system 
Collapse transformation matrix 

Loads .transformation that relates interface forces 
to the branch point load vector 
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Disc* fete deflections 
Discrete rotation 
Discrete accelerations 
Modal deflection 
Modal accelerations 
Applied forcing function 
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Interface acceleration function 
Shear stiffness 
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Abstract 

Thts docuraent details analytical procedures and dlsciisHOH 
example problem results for coupled base motion response analysis 
of payload structural systems. 

The report is divided into three major sections. 

I Analytical Development 
II Demonstration Example 
III Program Code 

A systems analysis program is described which by component 
analysis structural transient response analyses can be completed. 
Sophisticated aerospace structural systems require detail modeling 
and subsequent lot d and response analyses. The detail required 
consumes large amounts of computer time and a s' ^^nif leant numl>er 
of analysis samples. The subject analytical program presents a 
proven technique, used initially on the Sky lab program, which 
is designed to reduce cost and schedule time on detail structural 
analyses of Structural Payload Systems. 

Launch vehicles are composed of numerous structural components 
only one of which is normally defined as the payload structural 
system. However, to complete detail structural analyses, the 
entire structural model must be analyzed and the present custom 
is to reanalyze the entire system when the payload module or any 
component is revised, design model updated, or configuration 
updated. Detail or finite element modeling of large degree of 
freedom systems require significant amounts of computer time as 
well as complex models which are difficult at best for the engineer 
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or analyst to evaluate. Base motion procedures can be employed 
where critical segments of complex structural systems or com- 
ponents may be analyzed for various load conditions ithout 
having to re-establish the entire structural system coupled 
model properties. 

The transient response characteristics of a complex struc- 
tural system can be used as a basts for evaluating the transient 
response of a similar system, identical to the original except 
for changes in the dynamic representation of one or more of the 
original subsystems. 

The study is mathematically derived in Appendix A and 
illustrated step by step by a demonstration problem followed 
by a more detailed model typifying the HEAO-MECO event. 


f 
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I. Analytical Development 

A. Discussion 

Large complex structural systems create special problems 
for the structural dynamist. Launch systems, Atlas/Centaur, 

Titan, and Space Shuttle, have varied payload classes which can 
be accomodated as well as numerous mission events, conditions, 
and conf Igurs. t ions . These payloads require detail dynamics 
analysis to predict their design loads and accelerations. 

However, due to the inherent complexity of the structural 

systems, dynamic response analyses of the integrated system 

for each payload will necessarily be expensive and time consuming. 

The Skylab program had similar problems and requirements. 

Due in large part to the complexity of the docking latch and 
probe assembly a "base motion" method was derived to calculate 
orbital loads (Ret’s. 1, 2, and 3) using supplied docking tnt<.‘r- 
face forces which were implied to the latch connection points. 
Numerous analyses were completed improving the technique of 
applying interface or "base" forces to a structure culminating 
in Skylab Base Motion Analysis Report (Ref. 4), which defined 
loads analyses for Skylab launch configuration events. 

This culminated into the prime case motion technical approach 
(Refs. 5 and 6) by which incompatible base motion excitations 
are used. This led to the development, presented herein, of an 
analytical computer program for coupled bas j motion analysis of 
payload structural t^ystems. The methodology described In Appendix 
A was used to develop the associated computer program. Liberal 
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use of reference 6 was used to formulate this methodology. 


B. Technical Approach 

The primary function of the derived analytical program is 
to take the method described in detail in Appendix A, evaluate 
its application^ make modifications where required, and define 
a methodology whereby payload structural systems may be analyzed 
using a coupled base motion approach. 

This was accomplished by breaking the study into three 
distinct steps or cases for analysis. The subject cases compose 
three primary branch systems Identified in Figure 1. 


Main Support System, S Branch System, B Primed Branch System, B' 
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Figure 1 Primary Base Motion Systems 


Utilizing these' primary systems and applying the methodology 
contained within Appendix A a system of equations and techniques 
are derived which define an analytical method for coupled base 
motion response analyses of payload structural systems. 

Two cases are described to define a combination of payload 
effects and which make the basis for the derived methodology 
(see Figure 2). 
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Case A Evaluate response of the main support systems due 
to applied forcing function We now include 

structural effects of the branch system B and 
evaluate the total i-esponse from the known response 
of the support system, s. 

Case B Evaluate the response of the primed branch system, 
B' , by attaching it to the main support system, s. 
This is accomplished by removing the dynamic con- 
tribution of the branch system, B, from the total 
response . 

Case A is seen as the exact solution. Therefore Case B 


results correlated with the exact solution will yield the 
relative merit of the method. 
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The basic analytical approach of this study is to take 
the sample cases, Figure 2, apply these cases to a "hand crank" 
data evaluation using simplified data. Figure 3, to evaluate 
the analytical technique. The generation of the equations of 
motion for the subject study has been completed and is presented 
herein. A more sufficient demonstration model is included in 
Section II. 

CASE A 

Branch B representing the payload is defined by 
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The response equations are: 
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From the method suggested in reference 5, the equation of motion 
for this system subjected to base motion is: 

[“2b] {*i) ^ {*i1 + ['‘b] “ 

( 2 ) 

" [^2bJ 

and 

[%2j " [‘‘bJ ’ " L"’'bJ 

b] = - b2i] “ b] 

o°o {’'l) “ [®J K} + K} 
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Taking 

{*i}b - ^iJb Wb 


where 

[®iJb “ ^ 


and 

Jb " [V"2bJ " 

(3) 

or 

u) Jb - 94.8683 


then 

{^iIb - b] K} * [^j Wb 
K)b “ [®J + ['Plj Wb 

(4) 


and Mgg are subjected to inertia forces at the inter- 
face and to account for the total inertia force the two masses 
are lumped together as a part of the support system, booster. 


OKXGUNAL 
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where 


Mijj + Mgg “0,5 + 1.5-2 slu(5/ln 

- 2 slug/ia 

k « 1000 Ib/ln 
s 

For an undamped system without external excitation 
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From modal analysis we get: 



Thus Eq. (5) becomes 
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We also know 


( 8 ) 


{6<t)}={6j,g(t)]-+|6pj(t)j. and ■[S2(t)}=-[jc'2(t)]-+ii"2<t)} 


or 






where 


[“iJL - [0-5 O.5J 
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{*2<*>}“(*2<t>}+ Wg =rSE<*>}+ [“4, (4 


s s 


The interface forcing function (Fj^(t)} is determined by solving 
the booster dynamic equation. 
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Interlace Force-|Fj-(t^ 

For the branch, the response equations with interface 
force acting are 

(^1 “ ^2) " (^1 " ^2) “ “ 

( 11 ) 

“ZB^l * S C *1 ■ *2) + ‘'b (*1 “ *2) ■ “ 

• 4 

For inert! ally relieved system, Xg = 0, from (11) we obtain 

S (*1 - *2) + "b (’'l ■ *2) “ 

Taking [rj{Cg (ii - ig) + >^8 (*1 ‘ * 2 )} " 

and from Eq. (11) we get 

{Fj(t)}= - [r] [Mggj {ij 

Substituting for Eq. (4) obtain 

“ - [rj [mbb] ([»] (q}g) 

Substituting for Eq. (10) obtain 

= - [rj Kb] ([®J 

HEPEODUCIBILrrZ OP THE 
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and 

{Fj(t)}- [tb] {q} + [tsJ {^} + [tiJ 

where 

- - H [««] kl - - [>J [>] t‘] 

(1x1) ^ 

J - ■ H K.] C»J [-J - - [‘J [■] H [«■= 

(1x2) ^ 

- - [o.5 O.sJ 

.[-x] 


For payload or branch, from (2) using modal coordinates: 

- [“2b] [®] 

Substituting for Eq, (10) obtain 



( 12 ) 


(13a) 


0.5^ (13b) 


(13c) 


( 14 ) 
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and for the support or booster; 

{q}^ + {4}^ + {4^ - ['PsTlfift)} 


Substituting for Eq. (12) obtain 
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Eqs. (14) and (15) can be written in matrix form as 
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The various matrices are calculated below: 

Rj L'J [.,3. ■ [>] [•] [‘] [»■= «■=] 
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- [.;][ W [•] - - [■] [‘J W ■ t-‘J 



Substituting these into Eq. (16) finally gives 



The solution to this equation requires simply the definition of 
[6pjj(t)} which can be defined from the system described under 
1 of page 12. 
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Thus can be calculated solving Eqs, (17) and (18). 

The final deflection of the system can then be defined by 
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Finally {xj - [sj (xj + {<»}g 


- [ij K} [ij (ij (20) 
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where is found from Eq. (16) and '{xgj' from Eq. (19) 


CASE B 

Branch representing the new or revised payload. The 
base motion excitation is represented similar to Case A. 
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{xj - [bJ {X2} + {x=} 

Taking modal coordinates for Eq, (21) 

{xf; = 

Eq. (21) reduces to 

{q} + 2p«) j {4} + [" u)2j |qj. 
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['2p«J^,= [/kB,/M2g.j.[l.5492j 


and lumping the interface masses in the booster to define the 
support system we have as before: 


Total System 1 
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The response equation for Case 2 is; 
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M^g, + Mgg =2.0 slug/in. 
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kg = 10^ Ib/in. 

For undamped system without external excitation 
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From modal analysis; 
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{Fi(t)}> [tbJ + [t.sJ {5}^ + [tiJ {6pg} 


(29) 


Since M 2 g« * 1.5 Slug/in. we will have 
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(30b) 

- [0.75 0 - 75 J 
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Substituting these in Eq. (16) we obtain our final equation as 
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As previously noted in the response analysis of Case A (page 
IS) the solution to this equation requires simply the definition 
of [6pg(t)j. 
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From Case A we know: 
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To find 



the following steps are performed: 


(1) Analyze support of Case A with Fj(t) only acting 



(2) Subtract the results of Eq. (32) from the total 


response of support of Case A 
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be written as: 


where 
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Deflections obtained from step 2 for finding 

{«fe} 

Deflections obtained from solution of Eq, (31) 


and 

W - b] {*2} + ['i’JIb, W3. 

■ [^] {*2} + [o-sies] (qj 

- W + {<>i} 

where is found from Eq. (31) and -[xg} from Eq. (34). 


C. Conclusions 

During the course of evaluating the results, it was found 
that the details associated with modelling the interface required 


(34) 


(35) 
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further consideration. In particular, it is noted that applica- 
tion of the previously described approach converged to the exact 
solution when the mass on the branch side of the interface was 
made an integral part of the support system, s. That is to say, 
the branch side of the interface contains a stiffness element 
between the Interface and the first branch mass element. This 
fact does not impose serious restrictions on the basic approach, 
This effect is abridged by utilization of a massless interface 
for the branch side. 

The results obtained for Cases A and B are concisely 
summarized in Tables 1 and 2, and Figures 1 through 5. As 
indicated in the Figures for Case A, the curves are virtually 
the same. The same was also true for Case B, hence these 
figures are not included. Correlation of the two cases indicate 
errors less than 3%. Actually some individual errors were 
higher but seriously affected hy roundoff error. 

At this stage of development the method has been suffi- 
ciently validated to proceed with a more detailed analytical 
model which is defined in Section II. 
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3.091-1 


.159 

3.132-1 

.159 

3.130 


4.395-1 = 4,395 x 10"^ 


3 .352-3 
5.354-3 

2.986-2 

2.997-2 

7.674-2 

7.684-2 

1.080-1 

1.077-1 

1.342-1 

1.340-1 

1.641-1 

1.645-1 

1.840-1 

1.840-1 

2.137-1 

2.132-1 

2.619- 1 

2.620- 1 

2.922-1 

2.924-1 

3.052-1 

3.048-1 

3.346-1 

3.348-1 


.180 

3.348-1 

.180 

3.346 

N.B. 


3.934-2 
3 . 932-2 

7 050-2 
7.044-2 

7.453- 2 

7.454- 2 

7.961-2 

7.985-2 

1.134- 1 

1.135- 1 

1.688-1 

1.685-1 

2.178-1 

2,176-1 

2.329-1 

2.332-1 

2.309- 1 

2.310- 1 

2.591- 1 

2.592- 1 

3.190- 1 

3.191- 1 

3.661-1 

3.657-1 


.180 

3.348-1 

.180 

3.661-1 

Units: - 
Time: s 

Accelera 
Velocity 
Displace 



•’a 

1 ....... .1 

8.104-7 

1.923-6 

2.040-5 

2.038-5 

4.181-5 

4.494-5 

2.092-4 

2.095-4 

4.770-4 

4.786-4 

1.012-3 

1.014-3 

1.929-3 

1.927-3 

2.421- 3 

2.422- 3 

4.270-3 

4.269-3 

4.228-3 

4.225-3 

6.599-3 

6,598-3 

6.472-3 

6.471-3 

8.776-3 

8.776-3 

9.091-3 

9.096-3 

1.155- 2 

1.156- 2 

1.204-2 

1.204-2 

1.526- 2 

1.527- 2 

1.561-2 
1 . 560-2 

1.966-2 

1.965-2 

1.980-2 

1.980-2 

2.433-2 

2.433-2 

2.428-2 

2.428-2 

2.896-2 

2.896-2 

2,905-2 

2,905-2 

.180 

2.896-2 

.180 

2.905-2 

.180 

2.896-2 

.180 

2.905-2 


2,488-4 

2.487-4 

9.386-4 

9.380-4 

2.055-3 

2.054-3 


.180 

2.964-2 

.180 

2.966-2 
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TABLE II 
Case B 



(in) 

*2 

(in) 


*3 

(in) 


Exact 

Soln. 

Proposed 

Soln. 

1% Error | 

Exact 

Soln. 

Proposed 

Sola. 

1% Error | 

Exact 

Soln. 

Proposed 

Soln. 

1% Error ' 

j.015 

4.789x10“^ 

-1.523x10"^ 


8,532x10"® 

8.601x10"® 

■■ 

1.474x10"^ 

1.474x10"^ 


1 

j.030 

3 . 176x10"^ 

2,902x10“® 

9.44% 

2.084x10"^ 

2.086x10"'* 


9.377x10"'* 

9.386x10"^ 

0.1 % 

1.045 

3,506x10"*^ 

3.499x10"* 

0.20% 

-4 

9.937x10 ^ 

9.923x10"^ 

0.14% 

2.051x10"® 

2.053x10"® 

0.1 % 

.060 

1.507x10“^ 

1 . 501x10"® 

0.40% 

2.289x10"® 

2.287x10"® 

0.09% 

3.151x10"® 

3.152x10"® 

1.03% 

^0?5 

3.622x10“^ 

3.622x10"® 

0.0 % 

3.826x10"® 

3.829x10"® 

0.08% 

4.415x10"® 

4.412x10"® 

0 . 07% 

.090 

6.020x10"^ 

6.018x10"® 

0.03% 

5,808x10"® 

5.811x10"® 

0.05% 

6 . 147x10"® 

6.142x10"® 

0.08% 

i 

’,105 

8.223x10"^ 

8.227x10"® 

0.05% 

8.331x10"® 

8,326x10"® 

0.06% 

_3 

8.610x10 

8.611x10"® 

0.01% 

Li20 

1.059x10"^ 

1.059x10"® 

0.0 % 

l.llOxlO"® 

—2 

1.110x10 ^ 

0.0 % 

1 . 183x10"® 

1.183x10"® 

0,0 % 

4 

.135 

1.361x10"^ 

1.361x10"® 

0.0 % 

1.420x10"® 

1.421x10"® 

0.07% 

1.535x10"® 

1.535x10"® 

0.0 % 

.150 

1.742x10”^ 

1.741x10"® 

0.06% 

1.793x10"® 

1.793x10"® 

0.0 % 

1.878x10"® 

1.878x10"® 

0.0 % 

;.i65 

2.188x10"^ 

2 . 188x10"® 

0.0 % 

2.209x10"® 

2.208x10"® 

0.05% 

2.241x10"® 

2.241x10"® 

0,0 % 

1 

.180 

2.660x10”^ 

2.261xl0"‘^ 

'0 . 04% 

2,646x10"® 

2.646x10"® 

0.0 % 

_2 

2.676x10 

2.675x10"® 

1 0.04% 

.180 

2.660x10"^* 

\ 

0.04% 







i 

.180 


j 2.261x10"®^ 

1 


2.646x10"®^ 

2.646x10"®*' 

0.0 % 

— 2“* 

2.676x10 ^ 

-2* 

2.675x10 

0.04% 

1.180 












Average Error 

0.92% 

Average Error 

■■ 

Average Error 

0,04% 


* Maximum Values 


System Average 


0.34% 

















TABLE II (continued) 


Case B 


Time 

(sec) 

^1 

(in/sec) 


^2 

(in/sec) 


Xg 

(in/sec) 


Exact 

Soln. 

Proposed 

Soln. 

\% Error 1 

Exact 

Soln. 

Proposed 

Soln. 

1% Error 1 

Exact 

Soln. 

Proposed 

Soln. 

1% Error 

.015 

1.816x10”^ 

6.626x10“^ 


2.725x10"^ 

2.737x10"^ 


2.839x10"^ 

2.341x10“^ 


.030 

6 . 524x10"^ 

6.453x10"^ 

1.1 % 

2.971x10"^ 

2.968xl0"^ 

0.1 % 

7.040x10"^ 

7.050x10"^ 

0,14% 

.045 

4.317x10"^ 

—2 

4.348x10 

0.71% 

7.364x10"^ 

7.350x10“^ 

0.19% 

7.398x10"^ 

7.406x10"^ 

0,11% 

.060 

1.131x10“^ 

1.131x10”^ 

0.0 % 

9.486x10“^ 

9.503x10“^ 

0.18% 

7.534x10"^ 

7.510x10"^ 

0.32% 

.075 

1.599x10"^ 

1.596x10"^ 

0.19% 

-1 

1.134x10 ^ 

1.138x10”^ 

0.44% 

9.667x10"^ 

9.636x10"^ 

0.32% 

.090 

1.536x10"^ 

1.539x10"^ 

0.19% 

1.525x10"^ 

1.521x10"^ 

0.26% 

1.373x10"^ 

1.375x10"^ 

0,15% 

.105 

1.453x10”’^ 

1.456x10“^ 

0.21% 

1.792x10"^ 

1,788x10“^ 

0.22% 

1.918x10"^ 

1.921x10"^ 

0.16% 

.120 

1.762x10"^ 

1.756x10“^ 

0.34% 

1.915x10"^ 

1.921x10"^ 

0.31% 

2.313x10"^ 

2,312x10"^ 

0.04% 

.135 

2.281x10“^ 

2.276x10"^ 

0.22% 

2.264x10"^ 

2.266x10"^ 

0.09% 

2.326x10"^ 

2.326x10"^ 

0.04% 

.150 

2.782x10"^ 

2 . 789x10"^ 

0.25% 

2.674x10”^ 

2.666x10"’ 

0.3 % 

2.289x10"^ 

2.291x10"^ 

0.08% 

j.l65 

3.119x10"^ 

3 . 124x10“^ 

1.16% 

2.838x10“^ 

2.839x10"^ 

0.04% 

2.023x10"^ 

2.617x10"'^ 

0.23% 

1.180 

1 

3.123x10"^ 

3.121x10“^ 

0.06% 

3.032xl0”^ 

3 , 041x10"^ 

0.3 % 

3,176x10"^ 

3.170x10"^ 

0.19% 

j.l71 

1 

3.161x10"^* 

3.162x10"^* 

0.03% 











3.032x10"^* 

3.041x10"^* 

0.3 % 











3 . 176x10"^* 

3.170x10"^* 

0.19% 
_ 

1 

1 

1 

Average Error 

0.29% 

— 

0.22% 



i 

\ . . .. 

0.16% 


* llaxlmuB Values System Average = 0.22% 





















* Maximum Values 


*2 

(in/sec*) 



Proposed 

Soln. 

1% Error 1 

6 . 778x10"^ 

6.773x10"^ 


2.869 

2.856 

2.0 % 

2.371 

2.380 

0.4 % 

7.886x10"^ 

8,297x10^^ 

4.95% 

2.051 

2.019 

1.58% 

2.632 

2.573 

2.3 % 

8,941x10“^ 

9.564x10“^ 

6.51% 

1.322 

■ 

1.372 

3.64% 

3.074 

2.984 

3.01% 

_ 

1.928 

1.919 

0.47% 

6.896x10"^ 

7.919x10"^ 

12.92% 

2.143 

2.091 

2.5 % 


X 3 (in/sec^) 

Proposed 
Sola. 

3.447 
1.309 


2.940x10 >3.063x10' 
7.159x10"^ 6.931x10" 


3.325 

3.584 


2.102 

3.360 

3.563 

1.316 


7.242x10 *1-6.963x10 
8.021x10”^ 7.788x10" 
3.364 
3.565 


j% Error 

0.46% 
4.0 % 
3.3 % 

0.95% j 

( 

1 . 04% I 

0.59% 

1.52% 

4.0 % 

3.0 % 
1.46% 

I 1.85% 


^ a 

Ie 

td 

g*^ 

EO K 

m 


3.117* 


3.129* 


.4 % 


3.39% 


3.824=^ 


3 .820* 


System Average = 


0.1 % 
1 . 86 % 
2.62% 
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II. Demonstration Model 

A. Introductory Remarks 

A dynamic model is formulated to simulate a free~free 
vehicle in steady state accelerations. The transients of this 
model are examined due to application of external i) axial 
excitation and ii) axial and moment excitations. First, the 
complete modal system (i.e. integrated system) is analyzed as 
a basis for checking the results. Second, coupled base motion 
methods are employed to simulate the dynamics of the support 
plus branch combination to further evaluate the coupled base 
motion response analysis approach described in Appendix A. 

It cjhould be noted that for the latter we first require the 
interface acceleration response time history to augment the 
coupled base motion response equation, given as Eq. 14 in the 
aforementioned appendix. This can be conveniently obtained by 
dynamic response analysis of the support only for the external 
excitations acting on the support system with no branch, or by 
application of Eq. 17 of Appendix A if the previous system had 
an old branch. Comparisons of the results obtained for the 
complete modal system and fay the coupled base motion approach 
are made. 

B. Description of Model 

(1) Dynamic Model; Simplified structural models considered 
for the branch (or payload) and the support are shown in Figs. 9 
and 10 respectively. The weight, AE and El distribution assumed 
for these models are also shown in these figures. These 




3.0+10 1.5+7 22 1 Station 100 



SI Distribution AE Distribution Weight Distribution Branch Model 
lb“in^ lb. lb/ in 

TOTAL WEIGHT = 6600 lbs 

C.G. LOCATION: x CG (computed from z mass) = 109.1649 in. below Station 100 

z CG (computed from x mass) = 1.051 in. 


Fig. 9 Simplified Branch Model 


to 



El Distribution AE Distribution 
Ib-in^ lb. 



12.5 

Weight Support 
Distribution Model 
lb/ in 


TOTAL WEIGHT = 
C.G. LOCATION: 


5000 lbs 
X CG (computed 

z CG (computed 


from z mass) 
from X mass) = 


200 in. below 
Station 300 
0.0 in. 


Fig. 10 Simplified Support Model 
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distributions were basically obtained from similar data for the 
main axial beam of the HEAO-B Planar Model and Centaur Model. 
However, liberal simplifications were made with the aim of 
obtaining a simple demonstration problem. The total weight 
of the branch model is 6600 lbs and for the support it Is 5000 
lbs. An offset lumped mass of 3,63 slugs in the x and z plane 
i^ considered in the branch model, as shown in Fig. 9. This 
mass was transferred to joint 5 (Ref. Fig. 9) by the rigid arm 
branch coupling technique given in Appendix B. Thus only 9 
joints are considered in the free-free branch system. 

(2) Coupling Module: The Centaur Equipment Module is 

used to couple branch and support at the interface. The stiff- 
ness parameters of this module are given in Fig. 11. The coupling 
module was considered as an additional element of the cantilevered 
branch. Thus the cantilevered branch utilized in the base motion 
response analysis is as shown in Fig. 12. 


P, 6x 



Fig. 11 Coupling Module Stiffness Parameters 


46 



Fig, 12 Cantilevered Branch 


(3) Coordinate System and Degrees of Freedom (DOF): 

Figure 13 shows the coordinate system and the DOF considered at 
each joint of the models. 



y 


BBPEOBUCIBILrre OF TOE 

OMOWM. PAGE IS P(«l' 


Fig. 13 Coordinate System and DOF Orientation 


4 


Tables III and IV contain the coordinates and the DOF 
numbering system of the cantilevered branch and support joints 
respectively. 

TABLE III Cantilevered Branch Model Geometry 


Joint 

Number 

Coordinates 

DQF Numbering 

System 

X 

BH 


X 

z 

0 

y 

1 

100.0 

0.0 

m 

1 

2 

3 

2 

125.0 

0.0 

iil 

4 

5 

6 

3 

150.0 

0.0 

0.0 

7 

8 

9 

4 

175.0 

0.0 

0.0 

10 

11 

12 

5 

200.0 

0.0 

o 

• 

o 

13 

14 

15 

6 

225.0 

o 

• 

o 

0.0 

16 

17 

18 

7 

250.0 

0.0 

0.0 

19 

20 

21 

8 

275.0 

0.0 

0.0 

22 

23 

24 

9 

300.0 

0.0 

0.0 

25 

26 

27 

10 

300.0 

0.0 

0.0 

28 

29 

30 


TABLE IV Support Model Geometry 


Joint 

Number 


Coordinates 


DOP Numbering System 
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(4) Vibration Analysis; FORMA Subroutine VIB3 was used 
to formulate the stiffness and mass matrices of the cantilevered 
branch, free-free support and integrated system. This subroutine 
also gives the mode shapes, frequencies, loads and modal coupling 
transformations and collapse matrix (i.e, [ 3] ) . Table V contains 
the modal frequencies analyzed for cantilevered branch, free-free 
support and Integrated system. 


TABLE V Vibration Analysis Frequencies* 


Frequency in Hertz 

Mode No. 

Cantilevered 

Branch Free-Free Support 

Integrated System 

1 

5.84 

0 

0 

2 

28.08 

0 

0 

3 

54.31 

0 

0 

4 

79.40 

28.14 

7.78 

5 

126.68 

111.82 

27.85 

6 

153.15 

274.04 

44 . 7.3 

7 

250.86 

403.93 

80.13 

8 

314.23 

608.17 

99.28 

9 

337.72 

970.60 

115.50 

10 

475.06 

1511.70 

153.53 


*N,B. Data for 10 modes only given here. 


It was aimed to limit the frequencies to 100 Hertz so only 
4 modes of the cant ilevered branch system and 4 modes of the 
free-free support system were considered in the base motion 
response analysis. For the integrated system 8 modes were con- 
sidered. These mode shapes are plotted in Figs. 14, 15 and 16. 

Analysis for 8 cantilevered branch modes and 8 free-free 
support modes was also considered as a separate case. 
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Fig, 14 Cantilevered Branch Modes 


SEP 50 DUCIBIL!K OF THE 
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Freq. = 0 Freq. = 28.14 Hertz 

(c) Mode 3 Cd) Mode 4 

(Rigid Body) 


Fig. 15 (concluded) Support Modes Free -Free 
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Fig. 16 
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Fig. 16 (continued) Integrated System Modes 
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(5) External Base Excitation: Two cases are considered. 

Case A Axial excitation acting at joint 9 of support and 

along the negative direction of the 25th DOF. The 
variation of this base excitation with time is shown 
in Pig. 17. 

Case B Axial and moment excitations acting at joint 9 of 
support and along the negative directions of 25th 
and 27th DOP respectively. The variations of these 
base excitations with time are shown in Fig. 18. 

(6) Interface Response History: The support system was 

analyzed for both Case A and Case B excitations. The response 
equation used is; 

- fq3]g£p(t)} 

where Modal accelerations, velocities 

and displacements vectors for 
support 

2 

r<JuOg = Circular frequency squared 
matrix for support 

[<pjg = Support mode shapes (free-free) 
Pg = Damping coefficient for support 

Using the data obtained from FORMA Subroutine VIB3 for 
[cp]g and of the support and taking the damping coefficient 

pg = 0.01, FORMA Subroutine TR3 was used to formulate discrete 
interface acceleration history. 5 modes of the support were 
considered in this analysis. 


TIME THRUST 


SEC. 

LB. 

0 

30900 

.01 

30800 

.02 

30450 

.03 

29900 

.04 

27600 

.054 

60 JO 

.06 

3400 

.07 

2100 

.08 

1240 

.09 

880 

.10 

780 

.12 

0 

.50 

0 




Total Moment, ten thousands cf pounds- inch (dotted curve) 

Total Thrust, thousands of pound 


CD 


TIME 

SEC. 


26 


THRUST 

LB. 



0 

30 , 000 . 

24 

.042 

30 , 000 . 


.058 

29 , 756 . 

22 : 

.068 

29 , 300 . 


.072 

28 , 867 . 

20 . 

,078 

23 , 854 . 


.085 

17 , 500 . 

18 . 

.086 

16 , 020 . 

( 

.087 

14 , 738 . 

1 fi ! * 092 

r 

, .099 

1 

10 , 632 . 

5 , 402 . 

14 1 -102 

3 , 680 . 


.106 

3 , 160 . 

12 

.110 

2 , 700 . 

1 

1 

1 . 122 

1 , 920 . 

10 

1 . 142 

1 , 063 . 

j . 150 

711 . 

8 

1 .180 

0 . 


.50 

0 . 


6 




MOMENT 
LB. IN. 

8,000 

8,000 

9,780 

14 . 000 . 
18 , 940 . 

111 , 720 . 

212,000. 

212 , 400 . 

212,000. 

156 , 000 . 

67 , 200 . 

36 , 400 . 

27 . 600 . 

22 . 000 , 

15 . 600 . 
14 , 760 . 
14 , 220 . 
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C. Results and Observations 

Feed-Back Base Drive Program (see Section III) was used 
to formulate and solve the coupled base motion response equation. 
Some of the results obtained for discrete accelerations of the 
branch for Cases A and B are shown In Figs. 19 through 23. All 
these curves are drawn for results obtained considering 4 modes 
of support and 4 modes of cantilevered branch (4+4). Results 
obtained for the integrated system, considering 8 modes are also 
indicated in these Figures. As can be seen from these figures 
the results for the 4+4 modes, for both cases, agree very closely 
with the integrated analysis results. Results were also checked 
for 8+8 nodes (i.e. 8 support modes + 8 branch modes) and it was 
found that these fall between the curves for 4+4 modes and 
integrated system. It should be noted that in the 4+4 cases, no 
axial elastic support system mode was present (i.e. rigid body 
for X only), thus all elastic interface motion (axially) resulted 
from the feed back coupling. The curves are drawn for some of 
the degrees of freedom only, in both cases. Similar agreement 
of the results was found for other degrees of freedom. In con- 
clusion, we note that the method gave excellent results for a 
rather detailed model even with considerable modal truncation. 
Comparison of results from base motion methods with those ob- 
tained from the total integrated system analysis constitute 
sufficient arguments as to the validity of the coupled base 
motion methodology. 










Proposed Methodology xl3 
Exact k13 


xl3 (in/sec'^) 





Tyl5(rad/sec^) 



- 3 . 5 . 


Fig. 23 Case B - Tyl5 vs 


- Proposed Methodology Tyl5 
e Exact Tyl5 





III. Pro($i‘uiu Cotlo 

A. Program Logic 

Feed-Back Base Drive Program or FBBD is a program meant 
to formulate the coefficient matrices of the coupled base motion 
response eq.ua tlon and calculate discrete accelerations and internal 
loads of the branch system. The programming of FBBD is based on 
a collection of FORTRAN matrix algebra routines called FORMA 
(Ref. 7), as developed by the Dynamics and Load Section of the 
Martin Marietta Corporation, Denver Division. The user of this 
program is assumed to have a working knowledge of FORMA and its 
associated terminology. 

The program FBBD is divided into four main blocks, viz. 

(1) Block I: Formulation of the coefficient matrices of the 

coupled base motion response equation. The 
coupled base motion response equation used is 
derived as Eq. (14) in Appendix A and reproduced 
below. 
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[mIb 

[cp']b 

Ms 
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£qp] , tqgl , U3I 

{q'g 1 , {<ig i , Ug 1 


[tb] 

[r] 


[TS] 

[Til 

and [*6pg(t)} 


- Mode shnpoa matrix: of rnntllovomd branch 
= Mass matrix of cantilevered branch 

= Support free-free interface modes matrix 
= Branch circular frequency diagonal matrix 
= Support circular frequency diagonal matrix 
= Branch damping coefficient matrix 
= Support damping coefficient matrix 
= Collapse transformation matrix 

= Modal accelerations, velocities and displacements 
of cantilevered branch 

■ Modal accelerations, velocities and displacements 
of support 

=” Ratio of viscous damping and critical viscous damping 
(i.e. C/C^) 

= The interface loads transformation matrix which 
generalizes branch force to the interface 

= [B]’' 

- -[r][MlB(p][ffi']g ’ -Ib1’’[m1bIb!|5'1s 

- -miMlBfBl - -[sI'^ImIbIbI 

Interface acceleration distribution 


In the computer program FBBD Eq. (36) is defined as 
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(37) 
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or [A](q] + [B][q] + [c][q] = [d] [6^g(t) } 


where 


A(l,l) = 

A(l,2) - [fPclB^Mlg[0McP']. 


S 


A(2,l) 


[cp']^fe]MM]g[cp^]3 


= (^A(l,2) / 

A(2,2) = [l] - [cp']^[TS] = [l] - [<^>']gl^-[e]'^[M]g[0][cp']g) 


fl] + [P']s[e]'^[M]g[P][q^']g 


B 


B(l,l) = T2PUUJ 
B(2,2) - T2puujg 
c(i,i) = 

C(2,2) - r^j^jg 
D(l> “ 

and D(2) = [cp*]glTl] - -( cp' Igf P]^’l M]g[ Pi 


(38) 


The matrices [tp^.]g, [<P']g, fM]g, 

[P] are first generated using FORMA subroutine VIB3 
and are input In program FBBD. Data for damping 
coefficients, i.e. TpJg and Tp-Jg a£*e also input in 
this section. Using these matrices, the coefficient 
matrices [a], [b] , [c] and [d] are formulated as 
per the matrix operations defined in Eq . (38). 


Block I I : Solution of coupled base motion response equation 
Having formulated the coeffLCient matrices in 
Block I, Eq. (37) is solved for modal accelerations 
using subroutine TRIFBD (i.e. Transient Response for 
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Feed-Back Base Drive). This subroutine is basically 
a response routine to solve a second order differ- 
ential equation by fourth order Runge-Kutta (Gill 
modification) numerical integration. The initial 
conditions for velocity and displacements have to 
be established in order to solve this differ- 
ential equation. In our case, initial velocities 
at each DOF can be taken as zero, but the initial 
elastic displacements have to be computed using 
the data for external forces at initial time. The 
procedure to do this is explained in the next 
paragraph. 

The modal displacement vector Lq ] = 17 ,^- -j can 
be repartitioned, by placing the displacements 
connected with rigid body modes on top partition 
and the remaining elastic displacements below it, 
i.e. 

tq) = {- -q“-} (39) 


But 


and 




^ o 
% 

4. 


(40) 


== [ 0 ] 


t =0 


In view of Eqs . (39) and (40), Eq. (37) can be 
rewritten in partition form for t = 0 , as 
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The first of these equations Is 

f Mil] {q^ It^O" ^t=0 

■ • ^t-0 

The second equation of (41) thus gives 

[Mgi] [qj^lf^o ^t=0 

In view of Eq. (42) this reduces to 

^^e^t=--0 " ^^22^ (f^E^ ^t=0 “ ^^R^ ^ ^FE ^t=o) 

"" ^^ 22 ^ 1-^11^ ^^rO^^FE^^^ ^t=0 

(43) 

So in view of Eqs. (42) and (43) we can write 


1 

(44) I 







Eq. (44) can thus be used to find the displacements 


(>9 


at t “0, i.e. t!ie bottom partition oJ' tho de- 
sired initial elastic modal deflections. 

'\ 

0 


-- 


t=0 


(45) 


J t=0 


In subroutine TRlIilD the initial displacements 
are computed using the aforementioned approach. 

(3) Block III ; Calculation of discrete accelerations of branch. 

The discrete acceleration of the branch can be 
found using Eq. (2) of Appendix A, I.e, 


r6g(t)l - fPH6^(t)l + 




where {6(t)g} = Discrete acceleration of support 

= [cp* ]gt"6pg(t) ) + [cp']g[qg} 


(46) 


Substituting this in Eq. (46) yields 

0 


163 ( 1 )} 


\ P] [«P' ]gf^Fi;(t) } + 


[”c]b 


I [BHtp'l 





^b 1 


- ^ 

- 

Ri 


(47) 


from 


Having found the modal accelerations 
the solution of the base motion response “equa t ion 
in Block II and knowing the interface acceleration 

m ** 

history [6j,^(t)}, Eq. (47) is used to find dis- 
crete accelerations [6„(t>] of the branch, Thi.s 

c5 

is done in subroutine TACCBD (i.e. Total Acc(' Jo ra- 


tion of Base Drive). 
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(4) Block IV : Calculation of internal loads. The Internal loads 
ot the branch can be expressed ns 

{l^ 3 - |un'l|-Mjj^[6j3l . (48) 

Where ^1^3 ^ = Internal loads vector of branch 

[LCT] = Loads coefficient transformation matrix 
[M]g = System mass matrix of branch 

and = Discrete accelerations of branch, obtained 

from Block III 

The loads coefficient transformation matrix LlCt] 
and mass matrix Ls/llg can be obtained from FORMA 
subroutine VIB3 or by any other conventional struc- 
tural dynamics modal analysis package. In general 
a separate LlctJ for axial beam, bending beam (both 
shear and moment) and any spring stiffness elements 
are identified. The product [LCT][-M]g, thus gives 
the load transformation matrix which is discrete 
acceleration type. Based on this approach sub- 
routine ACCLD calculates the internal loads in 


br anc h . 
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B, DATA stream KlOW 


0000*000 «00*#00O000«00*«00*«0000O000*000<>*<HH*»*M**#'**tt«#000000*»«000 

DATA K1 / KB / K3 / 

0 70 / 15 / 40 / 

5 call start 

read ( 5*l00i) NCPB*NR0*NC0*Nul»NU2*NRTAPE*NXTAPEfNPRTfNWT*IFMC 

FORMAT (Ul5) 

R£aD{ 5*1002» IFINiTtTAPElD FORMAT <12A6) 

READt5»1003) sTARTTtOELTAT»ENDT FORMAT (10X*3E17.8) 

ReAO(5«1004) iFACCfIFLDS FORMAT (215) 

IF UFINIT .EQ. 6HINITIL) CaLL INTAPE (NWIiTAPEID) 

BLOCK I 


< 1) 


READ 

(A,NRTC,NCTC*K1,K1) 

(2) 


read 

(B*NRTItNCTItKl9Kl) 

(3) 


read 

{A*NRAfNCA»Kl»Kl) 

(4) 


READ 

{B*NRP8,NMPB*KltKl) 

(5) 


REaDIM (IVEC*NPI»NCI. ItKl) 

IF (iVECdJ ,NE. 999) GO TO 200 
DO 210 I - ItNRPR 


21 0 

IVEC(I) = 

I 

(6) 

200 

READIM 
GO To 110 

(DVEC»NRJtNCJ* 1 »K1) 

(7) 

100 

read 

(C*NCP8,NCTI.K1»K1) 

(8) 


READ 

(A»nCTI»NCTI*K1*K1) 

(9) 

110 

read 

(A»NDPO,NMPO»Kl »K1) 

(10) 


READIM 

tIVECfNRI»NCl*l*Kl) 

(11) 


reaOjm 

( JVECtNRj*NCJ,l ,K1) 

(12) 


read 

(2aR,NRZBtNCZB»ltKl ) 

(13) 


READ 

(QMBfNROB*NCOBf i*Kl) 

(14) 


READ 

(ZSUP»NRZS»NCZS»1»k1) 

(15) 


READ 

{0Ms*NROS*NC0Stl«Kl ) 

BLOCK 

II 


( 18 ) 


READ 

(TSEL*NRSfNCS»K2,K3> 


BLOCK in 


IF (IFACC .Eq, 0) GO To 500 

(17) READ (5*10fiii) ( STA ( I) » »NrP) FORMAT (12A6) 

BLOCK IV 


IF (IFLDS .to, 0) GO TO 500 

(18) REAUIM (LOs^NUNLDStlfKl) 

DO 495 L = liNLDS 

NLSS » LDS<L) 

(19) HEAD (C*NMOtMMOiKl#Kl) 

NL = 0 

DO 300 ILO = 1*Nl5S 
NL = LOS(ILO) 


REPRODlICii^- 
OtllGlKAL FAG 


UY OF TUE 
E iS POOR 


r>r»nr»f»noor» ortrt 


(20) REAO <B*NRB*NCB*KlfKl) 72 

(21) REaUIM (lVec*NRI*NCI*l*Kl) 

DO 290 I = l,NRI 

IF (IVeC(I) .NE* 0) ML * ML + I 
IF (NL .QT. Kl) GO TO 300 
290 CONTINUE 

(22) REaDIM tJVEc#NRj9NCJ*l#Kl) 

300 CONTINUE 

(23) READ (5*1002) ( STA ( I ) * I«i *N l ) FORMAT (12A6) 

a95 Continue 
500 continue 
GO To 5 
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C 

c INPUT matrices 
C 


ctl) 

A 

c 

COLLAPSE transformation 

C(2) 

a 

SB 

Interpolation transformation, branch restraint 

C 



Points to forcing function coordinates 

C(3) 

A 


Cantilevered branch mass matrix 

C{4» 

a 

■HP 

Branch cantilevered modes 

C{5) 

IVEC 

:s 

Identifies rows of preceeoing b to be kept 

C(6J 

JVEC 

M 

Identifies cols of preceeoing b to be kept 

C{7) 

c 

S 

Top partition (Negative) of d matrix 

C(8} 

A 

= 

Triple product of mass matrix and collapse 

C 



Transformation 

C(9J 

A 


Support free-free interface modes 

C(10) 

IVEC 

- 

Identifies rows of prfceeding a to be kept 

C(ll) 

JVeC 

*» 

Identifies cols of preceeoing a to be kept 

cnz) 

Z8R 

S 

bRANCH MODAL DAMPING VECTOR 

C<13) 

OMB 

s 

hranch modal Frequency vector (omega squared) 

C(14> 

ZSUP 

= 

support modal damping vector 

C{15) 

OMS 


SUPPORT modal frequency vector (OMEGA SQUARED) 

C(lfe) 

TSEL 

s 

tORCF selector MATRIX 

C(17) 

STa 

fe«t 

Response row identifier 

C(ia) 

LDS 

s 

Integer vector of size nlos 

c 



Elements s no. of sub loads to combine into 

c 



ONE load set 

C(191 

C 

= 

Mass matrix of cantilevered branch 

C»20) 

B 

= 

load coupling transformation 

C(21) 

IVEC 

s 

Identifies rows of preceeoing b to be kept 

C(223 

JVEC 

S 

Identifies cols of preceeoing b to be kept 

C(23l 

STA 

w 

Loads row identifier 




r'^ 


1 '.'* 


7 OF TiiE 
FOOR 
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C. Program Listing 

PROGRAM rtJBD (INPUTfOUTPUT*TAPES=INPUT*TAPE6=OUTPUT*TAPElOf 

* TAPEU*TAPEl2oTAPEltTApE2tTAPE7) 

DIMENSION A(70*70)* B{70f70)* C<70*70)» D{70#15)# 

• XVEC(70lt JVEC<70)» ZPR(70), 0MB(70)f ZSUP(70»* 
STA<70), TsEUa5*40)« VINIT(70J* DINIT<70)» LDSt70) 

DATA Kl / K2 / K3 / 

» 70 / l5 / 40 / 

DATA OMEPS / l-OE-03 / 


program to generate acceleration^ vELOCITVt AND DISPLACEMENT 
COEFFICIENTS FOR FtED-8ACK BASE ORiVE SYSTEM, 

solve for discrete ACCELERATIONS AND INTERNAL LOADS OF BRANCH. 


COUPLED BASE MOTION RESPONSE EQUATION* 
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DEFINE INPUT 

variables AND MATRICES 



NCPB = NO* 
NRO = NO- 
NCO - NO. 

OF BRANCH MODES TO 
OF SUPPORT I/F OOFS 
OF SUPPORT MOOES To 

BE KEPT 
(MAX a IS) 
BE kept 



NUl = SCRATCH tape UNIT 

NU2 s ScRaTcH TaPe UNIT 

NRTaPE = UNIT CONTAINING FORCING FUNCTION TJME HISTORIES 
NXTaPE - UNIT ON WHIcH OUTPUT HtSTORY WILL BE WRITTEN IN 

subroutine trifpd 
nprt = PRINT Interval for output 
NWT = forma Tape for output mATRICFS 
IFMC - 0* DISCRETE CANTILEVERED BRANCH 

= i» modally coupled cantilevered branch 

IFACC s If DISCRETE ACCELERATIONS OF BRANCH WANTED 
Of IF ABOVE IS NOT TO rE CALCULATED 
IFLOS = If inertial loads WANTeD 

Of If above is not to re calculated 


i 


r>ooonoor><ioo n ooonooooooooo or> oo 


START! = Initial time for output history 
END! = LAST tame FqR OUTPUT HtSTORY 
DELTA! 3 FORClNti FUNCTION TIME INTERVAL 

2RR = bRANCrt MODAL DAMPING VECTOR 

0M8 s branch modal frequency VECTOR tOMEGA SQUARED) 

ZSUP a SUPPORT MODAL DAMPING VECTOR 

OMS = SUPPORT modal frequency VECTOR (OMEGA SQUARED) 

TSEL a FORCE SELLEcTOR MATRIX 

DINIT a iNlTluL DISPLACEMENT (INITIAL ELASTIC DISPLACEMENTS 
ARE computed in SUBROUTINE TRIFBD) 

viNiT = initial velocity 

NLDS a NO. OF total load sets to CYCLE THRU 
LDS a integer vector OF SIZE NLDS 

elements a NO. OF SUB LOADS TO COMBINE 
INTO ONE load SET 

tool FORMAT (1615) 

1002 Format (i2A6) 

003 FORMAT tlOX,4El7.8) 

004 format (215) 

5 CALL start 

Read (5*1001) NCPetNR0*NC0»NUlfNU2rNRTAPE*NXTAPE»NPRTfNWTtIFMC 
READ (5*1002) l|: INIT.TAPEIO 
READ (5*1003) STaRTT.DELTAT.ENDT 
READ (5*1004) 1FaCC» IFLDS 

IF (IFINIT *£Q, 6HINITIL) CALL INTaPE (NWT.TaPEID) 


BLOCK I - FORMULATION OF COEFFICIENT MATRICES OF 


coupled oase motion response equation. 


read in branch mass matrix* collapse transformation* 
AND cantilevered MODES. 


c 


c 


C 

C 


c 


CALL HEAD (A*NRTC*NCTC*K1.Ki) 

call read (»9fNRTl.NCTI*Kl,Kl) 

CALL MULTB (A,B«NRTC*NRTI*NCTI*KI,k1) 


CALL write (B.nRTc.NCTI *6HT*MASS*Ki ) 

CALL WT^PE (a.NRTC.NCTI*6HT*MASS»Kl fNWT) 


CALL read (A»NRA»NCA*K1*KD 

CALL MULT (A*B9CfNRAsNCA»NCTl*Kl*Kl) 

CALL BTABA (A,B*NrTC#NCTI,k1,K1) 

WRITE (NxTAPE) i(A(l,j) *iai*NCTl) f J=l*NCTIj 


REPRODUCIBIUIT OF THE 
ORIGINAL PACK 13 POOR 


CALL READ IB.NRPe.NMPSiKlfKl) 

call zero (A.NRPb.NCPBiKI) 

CALL READIM (IVEC*WRI*NCI»IfKl ) 

IF (IVEC(I) .N£. 99«J) GO TO 200 
00 210 I a l*NRPp 
210 iVFCd) a I 

200 CALL READIM (oyEC f NRJ*NCJf 1 *Kl ) 

CALL REVaDO (l.tBf IVEC*JVEC*A,NRPB*NMPB*NRPBtNCP0.Kl*Kl) 


call WRITE 


(A»NRPOtNCPB*6HREnCM tKl ) 
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CALL WTAPL (A,NRPB»NCPB,6HREnCM fKltWWT) 

C 

CALL TRAnS !A,B*NRPBtNCP8fKltKl) 

Call MuLTB (ritC*NCPB»NRAtNCTl»KliKl) 

C 

CALL write (CyNCPBtNCTlvbH-DU) qKI) 

•C 

GO TO llo 

100 CALL read (C,NCP8tNCTI*KlfKi ) 

call read (A,NCTl,NCTIiKlfKi ) 

Rewind nxTape 
c 
c 

WHITE fNxTAPE) < (A ( I • J) * 1=1 iNCTl ) * J=1 1 NCTI ) 

C 

c read in support interface modes* 
c 

no call read ^a,ndpo»nmpo*ki*ki) 

CALL READIM (IVECfNRIiNCI»l*Kn 

CALL HEAOIM t JVECtNRjf NCJf 1 »K1 ) 

CALL ZERO (t>fNRO»NCO,Kl ) 

call REVADD (1.*A,ivEC,JvEC*B,NDP0*NMP0»NR0pNC0*K1vK1) 

C 

call write (tJ*NR0*NC0^6HRED0M tKl) 

CALL WTAPE £B*NRO»NCOf6HREDOM iK1,nWT3 

rewind NUl 

write (NUl) t ,I = 1 >NCPB) vJ=1*NCT1) 

, • C 

CALL MuLTA £C*8*NcPB»NR0tNC0,KliKl ) 

C 

CALL WRITE tC,NCP0iNCOt6HA(l,2» tKl) 

REWIND NU2 

WRITE (NU2) UC<I*J) I IsIfNCPB) tJ=lfNCO) 

C 

CALL TRANS (B,C*NROtNCO,Kl t Kl ) 

REWIND NXTAPE 

READ (NXTAPE) ( (A ( I #J) t I=I «NCTI ) » J=1 tNCTI ) 

REWIND NXTAPE 

CALL MULTB (Cf AtNcO»NCTl*NCTl*Kl*Kl ) 

CALL write (AtNCOfNCTI#6H-D(?5 »KI) 

C 

CALL MULTb (AiB*NcO*NRO»NCO«K1 *K1 ) 

call write (o*NC0»NC0*6hA ( 2t?) ♦Kl) 

c 

WRITE (NUl) ( (A(i»J) *I= lfNCO) #J=ltNCTt) 

WRITE <NU2) ( (B(I*J) ,I«1*NC0) .J=1*NC0> 

C 

C REaO in modal damping and frequency vectors# branch and 
c support. 
c 

call read (2BRfNRZa*NCZBil»Kl) 

Call read (dm0#nrob#ncobi1*k1) 

CALL read (2SUP#NRZS#NC2S#1,K1) 

CALL read <0mS#NR0S#NCOS.1#k1) 

c 

C FORM A» 8» C AND 0 MATRICES FOR EQ, OF MOTION. 

c 

Nmd = NCPb + NCO 

CALL ZERO (A,nMD»NM0#K1) 

CALL ZERO (B*NM0#NMD*K1 ) 

CALL ZERO IC,nM0*MMD*K1 ) 



o r» o on o o o n r> r> r> o o non 


77 


CALL ZEHO <0*NMD*NcTI*K1) 

DO 10 I = 1*NCP8 
A(I»I) a 1.0 

IF fOMB(I) .LE* OMEPS) OMB(I) = 0.0 
Bllyl) = 2.0 * * SORKoMpd)) 

Ca?!} = OMB(l) 

10 continue 

USE MEQ. FOR SUPPORT * 1.0 
DO 20 I a ItNCO 

IF (OMS(I) .LEo OMEPS) OMS(I) a 0.0 
B(I*NCPB.I*NCPti) = 2,0 * ZSUP(l) * SQRT(0MS<1)) 
C(I«NCPtif I+NCPBj s 0M5(Z) 

20 continue 


REWIND NUl 

READ (NUl) .islsNCPB) *Jes1*NCTi) 

READ {NUU ( (DU*nCPB.J) isleNCO) *J=lfNCTI> 

DO 30 I = It NMO 
DO 40 J = ItNCTi 
D(IfJ) = -D(ItJ) 

40 Continue 

30 CONTINUE 

REWIND IMU2 _ REPKoDUCIBILriY OF THE 

READ (NU 2 ) ((A(I,J*NcPB>*Ial*NCPBJ *J=ltNC0» ORIGINAL PAGE iS POOR 
Do SO I = 1 * NCPts 
DO 60 J a IfNCO 
A(J*NCPb,I) = A(ltJ+NCPB) 

60 Continue 
50 CONTINUE 

READ (NU?) i (A(l4NCPp*J+NCPB> tIal.NcO) »J«1»NC0) 

DO 70 I B It NCO 

A (I+NCPritl+NCPB) s 1.0 ♦ An+NCPBtl+NcPS) 

7o CONTINUE 


CALL WRITE 
CALL WTAPE 

call write 

CALL write 
CALL WTAPE 
CALL WRITE 


(A.NMDtNM0t6HC0EF,AtKl) 

CAtNMDtNMOtgHCOEFtAtKltNWT) 

<0tNMO»NMDt6HCoEF,a.Kl) 

(C,NMD.NMD»6HC0EF.CtKl} 

(C*NM0tNMOt6HCOEF,CtKl,NWT) 

iytNMDtNCTIt6HC0EF*DiKl) 


BLOCK i: - SOLUTION OF COUPLED BASE MOTION RESPONSE EQUATION. 


read selector matrix FOR DRIVING ACCELERATIONS. 

ZERO INITIL OISPLACEMENTSt READ IN INITIL VELOCITIES. 


CALL read (TSEL*NRStNCStK2tK3) 

CALL ZERO (OiNiT.ltNMOtl) 

call zero (VlNITtlfNMDtl) 


generate modal ACCELERATTONSt VELOCITIES. AND DISPLACEMENTS. 


REWIND NWT 

CALL TRlFBO (A.BfCtOtTSELtVINlT.OlNlTtSTARTTtOELTATtENOTt 
» KPRTtNM0tNRStNCStKltK2.K3tNXTAPE,NRTAPE) 


n o o o o o n o oo 
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BLOCK III - calculation OF DISCRETE ACCELERATIONS OF 3RANCH. 


IF (IFaCC *E0. 
REWIND NWT 
CALL ZERO 
CALL RTAPE 
call RTAPE 
CALL RTAPE 
READ <5tl002I 
CALL TaCCBD 


0) 60 TO 500 
CA»K1«K1«KI) 

(IARUN0«6HRE0CM t A*NRPt NCP»K1 «K1 tNWT) 
(lARUN0t6HRED0M t8«NR0«NC0»Kl tKl tNWT) 
<>ARUN0»6HT*MASS*C»NRRtNCR*Kl*Kl»NWT) 
(STA(I) *Ial»NRP) 

(A,B«C«STAf l.OiNRPvNRSvNMDvNXTAPEtNUU 
NPRT» STAPTT»ENDT*K1) 


BLOCK IV - calculation OF INTERNAL LOADS. 


IF <IFlDS .£«• 0) GO TO 500 
CALL READIM tLDS*Nl*Nl^'iSf l*Kl) 

00 495 L = 1*NL0S 
NLSS = LDS(L) 

CALL read <CtNMO,N«D*KltKl) 

call ZERo (A,K1*K19K1) 

NL = 0 

DO 300 ILD = 1* NLSS 
NL = LDSdLDJ 

CALL READ (tJ,NR8*NC8fKl fKl) 

call READIM <IVEC*NRI»NCItl<KD 

DO 290 I =! 1»NR1 
IF (IVEC(I) .NEo 0) NL = NL ♦ 1 
IF (NL -GT. Ki) GO To 300 
290 continue 

CALL REAOIM (UVEC*NRJ*NCJt 1 *Kl ) 

CALL REVADD ( 1 ,0 t B, T VECt JVECi A f NRB* nCB tNL*NM0*Kl iKl ) 
300 Continue 

CALL MULTA U«C*NL* W*NMD*Ki »K1) 

Oo 400 I = 1*NL 
DO 400 J ~ IfNNO 
A^fJ) = -AtI*U) 

400 CjNTINUF 
Pr'^^^ND NUl 

(5’1002) <STa<I» «I = 1*NL) 

call ACClD (AfSTA,sTARTTtENDT*NL*NMD»NPRT»KlfNUl) 

495 continue 
500 CONTINUE 
call LTAPE 
GO TO 5 
END 


<NWT) 
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subroutine TRlFao TSEL*XDO»Xn*STARTT*DELTAT»ENDT* 

NWRITEfNX*NF»NCS*KAfKF*KTfNTAPEfNRTAPE) 
dimension AtKAfil «B(KA*1> tC(KAtl) t D (KA« 1 ) t TSEL (KF • 1 ) * 

AC (42), FI { 15 J ♦F2tl5)»F(l5),oO(l30)fQ(I30j fP(4)t 
FC(7o), 

XDOa),XO(l) »XDO(130) •Xd( 130) tX(130) t 
XODMAX{130)» XDMAX(l30)* xMAX(130)< 

XODMir'dSO) • XDMIN(l30)t xMIN(130)» 
TDDMAXU30). TXDMAX(l30)» TxMAXa30)* 
TXDMINdSO). TXMIN(130) 


RESPONSE routine TO SOLVE 
(A)XOn * (8)X0 + (C)X 
FOURTH ORDER HUNGE-KUTTa 
IS USED, 

matrices a»b*c*u should not Share same core location (due 

THE ANSWERS (T«F|XDD»XD, X) WILL BE WRITTEN ON NTAPE EVERY 
ON PAPER EVERY NWRlTE * OELTAT. 

calls FORMA SUBROUTINES INVl *MULT»MULTB,PAGEH0« ZZBOMB, 
the maximum sizes ARE (BASED ON DIMENSIONS OF XODf XD*XtFtMAX,MINS) 
NX = 130 
NF = 15 


A 

B 


C 

D 


TDDMIN(130) « 
data NTT,N0T/5*6/ 

DATA NLPP / GO / 

C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
c 

G 

c 
c 
c* 
c 

2001 FORMAT «////15X»40H THE 

1 

? //23X« 


THE second order DIFFERENTIAL EQUATION 
= (D)F FOR XDO* XD» X, 

(GILL MODIFICATION) NUMERICAL INTEGRATION 


To MULTB). 
deltat and 


TSEL s 
XDo - 
Xo = 
STaRTTs 
DElTaT= 
ENDT s 

nwrite= 


(NXtNX) 
(NXtNX) 
(NXfNX) 
(NXtNF) 
I nCS ) • 


« 

« 


destroyed « 

DESTROYED « 
DESTROYED « 
DESTROYED « 


SUBROUTINE ARGUMENTS (ALL INPUT) 

= matrix coefficient of XDO. SIZE 

s matrix COErpiCIENT OF XO, SIZE 

= matrix coefficient of x. size 

= MATRIX coefficient OF F. SIZE 

FORCE selector MATRIX, SIZF (NF 

vectok of Initial velocities, size (nx), 

VECTOR OF initial DISPLACEMENTS. SIZE (NX), 

START time- 

integration STEP SIZE, 

END time, 

multiple of integration STEPS TO WRITE ON PAPER. 

NWHITE = 1 WRITE EVERY StEP (0,1*2#.,,) 

nwriTe = 2 Write every Sfcond step (0*2*4, •••) 

ETC 

SIZE of matrices A*e*C (SQUARE), NUMBER OF ROWS 
number of Cols in d. maxsIb 
number of cols in TSEL* MAX=40- 
ROW DIMENSION OF Af0*C*O IN CALLING PROGRAM, 

ROW DIMENSION OF TSEL IN CALLING PROGRAM, 

COL dimension Of tsel in Calling program, 
number of Tape on which answers will be written, (e.g* 


NX = 

NF s 

NCS = 

KA - 

KF = 

KT = 

NTaPE = 

NRTAPEs 


IN D 


The 


number of Tape 

DATA (TO 


10 ) 


11 ), 


xdd 

xo 


OUTPUT 

time 

fohce obtained by linear 
acceleration, size (NX), 

VELOCITY. SIZE (NX). 

DISPLACEMENT, SIZE (NX). 


containing driving forces. (E.Q. 

PE WRITTEN ON Paper and tape) is 

interpolation on TaBF. size (NF) 


INPUT Scalars to subroutine tri are • 

IOh STArTT s f 10.6, 
lOH DELtAT = F10.6, 




nno o o o orto 





-30 


80 

3 //23Xi lOH ENDt = Fl0.6e 

4 //^3X* lOH NWRTTE a T5 i 

?0i0 FORMAT t// TIME OFt 28^*SH TIME OF 

1 /«#Xi*»H HUWt 6XtlOH MAX AcCELt AXiiOH MAX ACCElt 

2 lOXflOM MIN ACCELi 6X*IoH MIN ACCEL 

3 // (lOX* I3» F16.6? E18.8* Flfl»6t E18.8) ) 

P020 FORMAT (// 20Xt8H TIME OFi 28X>fiH TIME OF 

1 /9X*AH Row* TXtSH MAX VEL* 8X*8H MAX VEL* 

2 12X*6H MIN Vel, 8^*8H MIN Vel 

3 // flOX* 13* F18.6, E18.8, FlA.6* El8,8) ) 

2030 format {// 20X*^H TIME OF» 28X*fiH TIME OF 

1 /9X*4H How* 7X*9H MAX DiSP* 7X*9H MAX OXSP* 

2 10X*9H MIN DISP* 8X*9H MlN DISP 

3 // (lOX* 13* F16.6, E 18 . 8 , F18.6* El8,8)J 
2040 FORMAT <//9X»8H TIME = F10,6) 

2050 FORMAT (//9X*15H APPLIED FORCES / (lOX* SElSeS) ) 

2060 Format (// 9X,4n row, 6X*13H acceleration* 8X,9H velocity* 

1 10X,I3ri displacement // dOX, 13* 3E20«8>) 


C 


NERROR*! 


IF (NX *GT. 130 .OR- NF ,GT. 15 ) GO TO 999 


compute initil displacements# 

REWIND NTAPE 

WHITE (HTAPEJ ( (a(I* J> *I=1*NX) *J=1*NX) 
write (NtAPE) ( (C ( I * J) * 1 = 1 *NX) * jsi *NX) 

Rewind ntape 

REWIND NRTAPE 

READ (NRTAPE) T» (FC ( I ) * 1 = 1 *NcS) 

Fc = interface accelerations at T = startt. 
CALL WRITE (FC*i.NCS#2HFC*l) 


CALL MULTW {TSEL*FC*NF*NCS*l*KFfKA) 
CALL MULTH <D»FC,NX*NF* 1 *KA*KA) 

Dn 70 J = 1 * NX 

IF (C(JiJ) #NE# 0 #) GO TO 70 

DO 65 r = 1 * NX 

CfI#J) = A{I,J) 

65 Continue 

70 CONTINUE 

call IN \/1 (C 9 A*Nx*KA) 

CALL MuLT (A*FC*X 0 *NX*NX. 1 *KA,KA) 

Read (nTape) ( ,i=i,nX) ,j=i,nx:< 
READ (NTaPE) nC(I,J) ,I= 1 ,NX) *J= 1 *NX) 
Do 80 I = 1 * NX 

ACC(I) = ACCd) ♦ RBMD{I*J) ^ F(J) 

IF {C(I*I) .EO. 0 #) X 0 (I)= 0.0 

80 CONTINUE 
Rewind nrtape 


i>CTBOT)'ICETLmr OP THE 
PAGE 13 FOOa 


PRINT INPUT SCALARS, 


call pagehd 

WRITE (NoT# 2001 ) STarTT,DELTAT*EMOT*NwRITE 
rewind ntape 

WRITE (Ntape) ((b(I*j)* i=i*nx), j=i*nX) 
call IN \/1 (A. B, NX* KA) 

write (NTaPE) ((b(I*J)* I = 1 *NX)* J= 1 *N> 1 ) 
REWIND ntape 


B=B 

BsAI 


onoo-o'on orjo oon 
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READ (NTAPE) I=XtNX)* J=1?NX) 

READ (NTAPE) I = lfNX), J=lfNX) 

call MULTB (At Bt NX, NX, NX, KA, KA) 

CALL MULTB Ui C, NX* NX, NX, KA, KA) 

CALL MULTb <A, Dt NX? NX, NF, kA, KA) 

REWIND NTAPE 

CALL write (BvNX,nX,6HINVA«B,KA) 

CALL WRITE (C,NA,mX,6HINVA*C,KA) 

CALL WRlTt (D,NX,fgF,6HINVA«D,KA) 


BsB 

A=AI 

BsAIB 

CsAIC 

D=AID 


SfT initial VALUES, 

NSTEPS = 0 
T s STaRTT 
DO 3o 1*1, NX 
Qd ( I ) = 0,0 
Q (I) = 0.0 
XO(I) = XDO(I) 

Jn X (I) - XO (I) 

REWIND NPTAPE 

READ (NRTAPE) n, (ACd) tI = l,NCs) 

CALL write (AC,1,NCS,2HAC,1) 

call mult (ISEL, AC*F1,NF,NCS»1 ,KF,KT) 

CALL write IFI vl,NFr?HFl,l) 

Do 205 I = 1 *NF 
205 F (I) = Fl (I) 

DO 38 1=1 ?NX 
XDD(I) = 0.0 
Do 37 J=1,NF 

37 XODd) => XDDU) * D(I,J)*F(J) 

DO 38 J=1,NX 

38 XDOa) = XDD(I) - B(I,J)*XD(J) - cd,J)«X(J) 
DO AO t-l,NX 


XDDMAX(I) 

- 

XDOtl) 

AC'DMlNd) 


XDD(l) 

TOOMAX d) 

S 

startt 

TODMINdl 


staktt 

XOMAX d) 

= 

XD(i) 

XPMIN(I) 

= 

XDd) 

TXDMAX(I) 

= 

START! 

TXDMIN( I) 

s 

staktt 

XMAX(I) 

= 

xd) 

XMTNd) 

= 

X( I) 

TXmAX(T) 


startt 

TXMIN(I) 


staktt 


SET INTEGRATION CONSTANTS. 


Ptl) = -5 

P(2) si.- SQRT(.S) 
P(3) s 1. + SQRt(.5) 
P(A) = .5 

WRITE DATA AT START, 


NW NWRITE 

eo TO 3*0 

Integration loop. ij=i,half step), (j=2,half step again), 

(J=3,FULL STEP), <J=A,END OF STEP). 

gill factor = .5 


So o o o o o o o 


S2 

100 READ (NRTAPE) Tit (AC ( I) *1=1 *NCS) 

C.'^LL MULT tiSEL*AC*F2*NF*NC5*l*KF*KT) 

DO 150 J=i*4 
00 no I = 1*NX 
2 = Xo (I) * DElTaT 

20 = XODCn « oeltat 

IF (J ,EQt A> 60 rO lOS 
R * P( JJ « (2 -Q (in 

rd = P(U) ^ (ZD"*oDan 

60 TO 107 

105 R « (7 “ 

RD = (ZD - 2*»00(I) )/6. 

107 X (I) = X (I) ♦ R 
Xnd) = XD(I) + RD 
0 (I) = Q (I) * 3.*R - P(J)<*Z 

110 QD(I) = QD(I) * 3,<»RD - P(J)*ZD 

GO TO (2o0»140f220» lAO) ♦ J 
200 DO 210 K = 1»NF 
210 F(K) = IFKK) + F2(K)) / 2,0 
GO To 140 

220 DO 230 ^ ~ 1*N? 

230 F(K) = F2(K) 

NSTEPS = NSTEPS + 1 

T = sta«tt + Float (NSTEPS) <»deltaT 
140 Do 15q Is 1 >NX 
XdD(I) » 0.0 
DO 145 KslfNF 

145 XDOm = XDDd) * 0(I*K)«^F(K) 

DO 150 KcIiNX 

l5n XDD(I) = XDOai - 8(I*K)»XD(K) - cU*K)«X(K) 

Do 240 I = 1 *NF 
240 FI tl) = F2 ( I) 

GET MAXIMUMS AND MINIMUMS. 

DO 330 I=l*NX 

IF (XDD(I) .LE. XDDMAX(I)) GO To 30S 
XDOMAX(I) = XDJ(I) 

TDDMAX(l) = T 

305 IF (XDD(I) *6E* XDDMlNd)) GO TO 310 
XDDMIN(I) = XDO(n 
TDDMIN(I) = T 

310 IF (XDtn *LE. XdMAX(I)) GO TO 3lS 
XOMAX(I) = XOU) 

TXOMAX(I) = T 

315 IF (XO(I) .GE. XOMIN(I)) GO TO 320 
XOMIN(I) = XO(I) 

TXOMIN(I) = T 

320 IF tXlI) .LE. XMAX(Il) GO TO 325 
XMAX(I) = X(I) 

TXMAXd) = T 

325 IF (XI I) .GE. XMlNCin GO TO 330 
XMIN<T) = X!D 
TXM ■ = T 

33{j COf 

WRlTE ANSWERS ON NTaPE FOR SUBSEQUpNT USE (SUCH AS 

Time response addiIional equations or plot). 

340 write (NTAPEJ T# (F(J)* J=1«NF)* ( XOD ( i ^ * XD < I) *X ( I) # I=1*NX» 
SEF data should PF prtmtfd. 


r>r>o ooo non 



^ IF (NW .LT. NUf^iTEJ qO to 34S 

i CALL PAtiEHD 

WRITE (NOTf2040) T 

! write (NOT^iiOSO) (F(I), I=l,NF) 

1 WXS = 1 

jl . NXE S NX 

NFLN = tNF-U/5+1 

i IF <(NX£ * NFLN) ,GT. <NLPP-iS) J NXE= {NLPP-15) -NFLN 

342 WRITE (NOT*cObO) {I* XOO < I ) » XD(D* X(I)f I=NXStNXE) 
1: IF (NX .EQ. NXE) GO TO 343 

NXS = NXE * 1 

I NXE = NX 

IF ((NXE-NXS) eOT. (NLPP- 9)) NXEaNXS+ (NLPP- 9) 

I CALL PAGEHD 

GO TO 342 

I 343 NW s 0 

345 NW = NW+1 

SEE IF RUN HAS DIVERGED, 


NERR0R»2 

DO 350 I«1*NX 

IF {AB5(X(D) ,6T, l,E + 35) GO To 999 
350 continue 


DETERMINE IF RUN IS FINISHED. 
IF (T ,LT, END'D GO TO 100 


PRINT MAXIMUMS AND MINIMUMS IF RUN IS FINISHED* 


00 4i0 MM=1»3 
NXE = 0 

400 NXS = NXE + 1 
MXE = NX 

IF ((NXE-NXS) *6T. (NLPP-ID) NXE=NXS+ (NLPP-1 1 ) 


CALL PAf^EHD 

IF (MM *EU* 1) WRITE (NOT.2010) 

IF (MM *EQ. 2) write (NOT*2020) 

IF (MM *EQ* 3) write (NOT*2030) 

IF (NX *GT. NXE) GO TO 400 
410 CONTINUE 
RETURN 

999 CALL ZZ90MB (6HTRIFBD*NEPR0R) 

end 


(ItTDDMAX(I) fXDOMAX(I) * 

TDDMIN(I) «XDDMIN(I) • ZsNXStNXE) 
(IfTXDMAX(I) tXOMAX( I) « 

TXnMIN(I) *XDMIN(I) » I=NXS*NXE) 
(ItTXMAX(I) #XMAX(I) * 

TXMIN(I) *XMIN(I>» IsNXSfNXE) 


™ AL IS POOR 
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SUBROUTINE TaCCdp < PHlBR*PSUPfRRMDtSTft>CONVT»NRP*NF*NX* 

NRTAPe,NXTAPEflNPRlNTtSTARTTfENDT*KR) 
dimension PHZBR(KR*M« PSUP{KRM)f RRmD(KR»1)o STA(1)« 

* Fa5)» XDO{70)« ACC(7o)« ACCMAX {70*4} i TlME(S)* 

# THlSTSTOtS) 

C 


PROGRAM TO GET TOTAL ACCELERATIONS FoR BASE-0R|VEN PAYLOAD 
TOTACC = PHI (CANT) ^ BDACC 

* PHKrB) « PHI (SUPPORT SYS) * SUPPORT ACC 

♦ PHI(RB) BASE ACC 

DEFINITION OF INPUT VARIABLES 

NRP = NUMBER OF RESPONSES TO BE CALCULATED 

(Number of rows in phibr) 

NF - number of PSEUDO-FORCES USED TO DRIVE MODEL BASE 

NX = number of total MODES CONSIDERED IN THE RESPONSE EQUATION 

OF THE coupled SYSTEM 

NRTAPE = logical UNIT OF DRIVING RESPONSE TAPE 
NXTAPE = LOGICAL UNIT OF OUTPUT RESPONSE TAPE 

nprint = PRINT Interval for response calculations 
startt = start Time for response calculations 
ENOT “ END time FOR RESPONSE CALCULATIONS 
ONVT = scale Factor for answers 
PHIBR =* CANTIlLVEREd MODES FOR DRIVEN SYSTEM 
PSUP = reduced free-free modes of DRIVING SYSTEM 
RBMO = REDUCTION TRANSFORMATION FOR DRIVEN SYSTEM* 

relative to base 

STA a RESPONSE ROW IDENTIFIER 


1002 FORMAT (10X*4El/,8) 

1003 FORMAT (l2A6) 

aooi format c//56x,#acceleration Time history*//iox»<»sta ** 

» 5(7X»»TIME = *fF8.5)/) 

2002 Format (2X9l592X,A6,S(7XfEl5,8 ) j 

CALL MULT (R8M0*PSUP*PHIBR(i*NCP+i) fNRRiNRO«NCO*KR*KR) 

call write (Phibr, nrp*nx*6hc,qdd ,kR) 

REWIND NR?ApE 
REWIND NXTAPE 

read PREVIOUSLY WRITTEN TAPE AND CALCULATE ACCELERATIONS 

NTIME = 0 
iPRINr a 0 

UO READ (NRTAPE) T* (F ( J) # Jsl *NF) * (XDD ( I) ,XC » Xt I«1 *NX> 

CALL zero (Acc*NRP,1,KR) 

DO 20 I = 1*NRP 
DO 30 J s l#NF 

ACC(I) = ACCd) ♦ RBMD(I*J) ♦ "(U) 

30 continue 

DO 40 J = If NX 

ACC(T) ss ACC(I) + PHiBRdfJ) # xOD(J) 

40 Continue 

ACC (I) = ACC (I) * COnVT 
20 CONTINUE 


o o ri o o r> 
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TT •• . ,v,- 

“^•1 

i’-l 


IF (NXTaPE 0) WRITE (NXTaPe) Tt (aCC ( J) * J^i »NRP) 

SEE IF MAX OR «IN hAS SEEN CALCULATED AT THIS TIME 


IF (T •UT« STAKTtJ go TO 50 

CALL INTMAX " < ACCf ACCMAX*NRPfT,KR) 

GO TO 6U 

50 CALL maxi (ACC*ACCMAXiNRP*TtKR) 

60 CONTINUE 


PRINT DATA 


IPRINT = IPHINT * 1 

IF (IPRINT *£00 NPRINT .OR. T ,fQ, STaRTT .OR. T .EQ. ENDT) 
* GO TO ?o 

GO TO 80 

70 NTIME a NTIME * 1 
IPRINT = 0 
TIME(NTIME) = T 

CALL ASSEM (ACC*1,NTIME*THIST»NRP.19NRP»5*K1*KI) 

IF (NTIME •EO^ 5 .OR, T .EQ. ENnT) GO TO 90 
GO TO flO 

90 CALL PAGEHD 

WRITE (o,200l) tTIME(J) *JslfNTlME) 

NLINE r 0 
DO 100 I S 1*NRP 
NLINE a NLINE ♦ \ 

WRITE (P*2002) I,STAfI) ♦ (T hIST(i.J) fJalfNTIME) 

IF (NLINE .lT, Ao> go TO lOO 

CALL PaGEHO 

WRITE (fe.aooU ITIME(J) *JsifNTlME) 

NLINE a 0 

100 continue 

NTIME a 0 

80 IF (T ,LT. ENDT) GO TO 110 

CALL WRITE (ACCMAX*NRP*4»6HAcCmAX.kR) 

RETURN 

end 


I 
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SUBROUTINE INTMAx (A*B*NRA* VAR*kRA) 

DIMENSION A(l)»tf{KRAa) 

C 

C this program INlTIAtiZES MATRIX B WITH THE VALUES IN VECTOR A AND 
c The variable var 

INPUT VARIABLES 

A = VECTOR OF VALUES TO 8E STORED IN COLUMNS 2 AND A OF B 

8 s OUTPUT Resultant matrix 

NRA S number Cf elements IN A 

VAR 3 VARIABLE TO BE STORED IN COLUMNS 1 AND 3 OF B 

KRA = DIMENSION SIZE OF A AnD ROW SIZE OF B IN CALLING PROGRAM 

DO 100 IRA=1«NRA 
B(IRA*! )=VAR 
B(IRAf3J=VAR 
B(IRA»2)s;AaRA) 

BdRAfAi-sAIIRA) 

100 Continue 
Return 

END 







I 


oononooooooonooon 


SUBROUTINE MAXI ( A*BtNRA# VAR*KR aJ 
niMENSlON AtDfBcKRA,!) 


THIS SUBROUTINE COMPARES VALUES IN VECTOR A WITH THE VALUES 
In columns 2 AND 4 OF MATRIX B 
IF (A tl) .GT,B(lf2) ) Ba,lJ=VAR 

B(I,2)aAm 
IF (A(I) «LT*B(I*4n Ba*3)=VAR 

B(I«4)3A(I) 

INPUT variables 

A = VECTOR CONTAINING VALUES TO BE COMPARED WITH VALUES 
JN S 

B = (NRA»4) Matrix of max and mins 
NRA S number of ELEMENTS IN A 

VAR s A constant which REPLACES VALUES IN COLUMNS 1 AND 3 OF 6 

KRA a dimension size OF A AND RoW SIZE OF B IN CALLING PROGRAM 

DO 110 1RA*1»NRA 

IF{6(IF. f2),GT,A(iRA>) GO TO lOO 
B{ TRA,l JaVAR 
B (lRA*2)aA (IHA) 

100 if{b<ira»4) .lt.a(ira)) go TO no 

B(IpA,:i)=VAR 

0(TRA94)=A(IRA) 

110 continue 
RETURN 
END 


non ooo oono-^ooo o o o o o o o 
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SUBROUTINE ACCLO { A»STA*STARTT*ENDTt N(. *NX*NPRINT»KA«NRTAPE) 
DIMENSION A(KA,1) , STAID * ZLDiTo) •ZLD mAX ( 70»^) t 
* TIME(&), ThIST(70»5) 

Data kz / 70 / 

this subroutine Calculates internal loads in branch 


INPUT variables 

A = loads transformation (DISCRETE ACCELERATION TYPE) 

STA 3 load row identifier 

STARTT a start TIME OF RESPONSE CALCULATED FOR BRANCH 
ENDT = END time OF RESPONSE CALCULATED FOR BRANCH 
NL = NUMBER OF ROWS OF MATRIX A 

NX 5 NUMBER OF COLS OF MATRIX A 

KA = ROW OAMENSION OF A IN CALLING PROGRAM 

NRTapE = NUMBER OF TAPE CONTAINING DISCRETE ACCELERATION 
DATA OF branch 

NTIM = 0 
IPRINT = 0 

UO READ (nRTAPE) T» (2LD( J) *J=l*NX) 

CALL MULTB {A,2L0,NL*NX,1 *KA,KZ) 

SEE IF MAX. OR MIN. HAS BEE CALCULATED AT THIS TIME 

IF (T ,GT. SfARTTJ 60 TO 50 
CALL INTMAX ULD.ZLDMAX.NLfTtKZ) 

GO TO 60 

So CALL Ma^I <ZLDiZLDMAX,NL»T,kZ) 

60 CONTINUE 

Print data 


c 


c 


c 


IPRINT = IPRINT ♦ 1 

IF (IPRInT .EQ. NPRSNT .or. T .eQ. STaRTT .OR. T .EQ. ENOT) 
♦ GO TO 7o 

GO To 80 


70 NTIME = NTIME ♦ l 
IPRINT s! 0 
TiME(NTImE) = T 

CALL AsSEM (ZlO»1, NTIME. THIST*NL*1.NL»5*KZ*KZ) 

if (NTiME «E0. 5 ,OR. T ,EQ. ENdT) 60 TO 90 
GO TO 80 


:all pagehd 

miTE (6.2001) (TIME(J) .J=1 »NTImE) 
JLINE = 0 
)0 100 I a 1»NL 
JLINE = NLINE ♦ i 
(RITE (6,2002) I .STA ( I ) • (THIST ( I 
:f (nltnf .lT, AOS fln To ion 


T),mor>'’^CETI.ITY OF raE 
Snae page is poor 


I) .1.1 aMTTMPt 


call PAGEHD 

WRITE (6,2001) (TIME(J) .Jsl.NTiME) 
NLINE s 0 

100 continue 










80 IF «f .Lf. ENOTi 60 TO 110 

CALL write (^LDMAx*^^L»4,6HZL0MAX•KZ) 

RETURN 

END 


. .... , ■ 


APPENDIX A 

Salient Mathematical Features of the 
Coupled Base Motion Response Technique 


Saltont Mathematical Features of the Proposed Method : 

Schematically the basic method to determine the dynamic 
response of any branch system using the response analysis of the 
interface, obtained from the response analysis of an integrated 
system (base-branch) incorporating a different (or no) branch, 
can be depicted as shown in Fig. lA. 

For malhematical convenience, the derivations in part (b) 
of Fig. lA, namely the coupling of an alternative branch to form 
a new integrated system, will be presented first, followed by 
part (a). 

The transient response equation for the branch subsystem 
denoted herein as component B, subject to interface motion 
excitation only, can be written as the following two dynamical 
subsystem equations (reference 5); 

[M]g{x^}g + [C]g{x^}g + [K]glxj^}g 

( 1 ) 

=-[M]gf0][6(t)} 
and [x}g - [0](6(t)} + 

where " Displacement of branch B in complete system 

{x^]g = Displacement of branch B with constrained 


Interface displacement 



Fig. lA 


Schematic Diagram of the Proposed Method 
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{6(t)} - 


Interface displacements of the support system 


[ 3 J " Collapse transformation 

« Mass matrix of constrained branch 
B 

cl ° Partitioned damping matrix 
[K ” Partitioned stiffness matrix 


Physically, Eq. (1) represents the constrained branch, and Eq. (2) 
represents the free branch of the dynamical subsystem. Applying 
the transformation equation 


where 




qr - Vector of normal modal coordinates of con- 
strained branch 


(3) 



Cantilevered mode shapes, normalized to 
yield unity generalized masses 


and substituting Eq. (3) into Eqs. (1) and (2) yields 


{q} + ['Spuj J {q} + {q} 


- - ['f’clg Mg b] 


R K?K^ *3 ) I A ■ i . . • ; i i'Y ( 
ORIGINAL 13 POOR 


(4) 


■ • 


K Wfcai ft --' 
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Wg “ [s] {^<t>} + [^clg Wg 

where tu =» circular frequency 

p *= Ratio of viscous damping and critical viscous 
damping (i.e. 

Equations (4) and (5) would be directly applicable to the pro- 

•• 

blera only when the interface motion {6(t)) was derived from an 
integrated system analysis incorporating a branch model identi- 
cal to that used in Eq. (4), ' 

The new branch (or component B) when coupled to the base 
or support (referred to as component S) will result in new inter- 
face forces differing from those resulting from the original 
branch load-class. The interface forces are used to couple the 
new branch to the base system in the interface motion response 
equation. Finally the situation as shown in Fig.2A is obtained. 





{Qpg(t)]* 

(External forces act on 
Component S only) 


(FjCt)} 


Interface Forces 


Interface 

Accelerations 


Fig. 2A 


The base motion accelerations input to part (b) of Fig, lA 
cori'espond to the acceleration of the interface coordinates of 
component S resulting from the external forcing function {QggCt) 

f 1 

only; consequently we denote this interface motion as 



Similarly, we denote the acceleration of the interface of the 


coordinates of component S (base) under the action the intor- 
I 

response analysis into 2 parts as shown in Fig. 3 & 


f 

face forces *^F 


(t)! as Thus, we can divide the base 


Base System Part 1 . Part 2 

^ /\ ^ ^ A j ^ A 



Fig. 3 a 


The application of superposition principle then yields 

C6) 

At this stage, it is worthwhile mentioning that in defining 
the dynamic discrete model of the base system, interface masses 
of the branch will have to be lumped with the corresponding 
interface masses of the base. Thus the full inertia of the 
interface masses is accounted for. In the interest of separating 
branch analyses from support system analyses, it is suggested to 
use an interface branch point with no mass. 
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The response equations of the base system under the action of 
» 

interface forces {^jCt)} (Part 2 of Pig.3A) is shown in Eq, (7). 


{q} + {q> + J {q}. 


W! 


(7) 


where j^cpj * mode shapes of component S normalized to 

yield unity generalized mass matrix. 

Thus the interface acceleration due to (Fj(t)j‘ can be expressed 


as 
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{spiCt)} = [cr'l {5} 


(8) 


where 


“ The pertinent rows of matrix 

to the interface coordinates only. 


corresponding 


Substituting Eq, (8) into Eq. (6) yields 
{*«)} - + ['P’lg {«}s 

The interface forces {Pj(t)} can be computed from the inertia 
forces of the branch as follows ; 

{Fj(t)} = - [r] [m] {i} (10) 

B B 


whore j^pj - matrix that relates inertially relieved 

forces to interface forces. 

Substituting Eqs. (5) and (9) into (10) yields 

{-,<•>) ■ - [r] H, {[.] * [«] [•■], «), 

• W, W,) 

which can be rewritten as 

{Pj(t)} ” [tb] {q} + [ts] {q} + [ti] {iSjg} 

B “ 


OKUnS.'F- . 




rOOA 


( 11 ) 
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Where 

[tb] - - [r] [m]^ 

[ts] = - [r] [m]^ [ 9 ] [cp-]^ 


[-] - - [-] [«; [ 3 ] 


Substituting Eq. (11) into Eq. (7), yields the response of the 
base in the coupled system: 

{5} + { 4 } + {q} . [ctff ([tb] {4} 


B 


[^s] {4} + [ti] {6^3}) 


( 12 ) 


The response of the branch in the coupled system is obtained by 
substituting Eq. (9) into (4), yielding: 


{4}^ H. {4}^ + {q}^ - 

- [<^ 0 ]; W3 [ 3 ] {3 bs} 

- Ml [m; [3] [.•; {^}3 


Eqs. (12) and (13) represent a system of simultaneous - 
coupled equations which can be combined into one expression 


as follows : 
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] I 
1 

w: 

; [“]g [®] [ 

w’l 
■ -'s 

_ 

[•C ["] i f ■ j 
^ 1 

- b 

'3g C’f' 

*1 


r- ^ 

1 ~ 

> + 

2p>", 


0 


- 


• » 


0 

i: 




% 
- ^ 



r 

d 

1 r 

B 

H [a]' 


q 

— — >■ SB 

s 

u 

[-] 




(14) 


Thus the transient rei ponse of the new system with an 
incompatible interface motion input requires the solu- 

tion of Eq. (14) 

With respect to part (a) of Fig. lA the response of the 
original integrated base- branch system is known, i.e. we know 
■[6(t)j- and {^^(t)^'. The response of the base system with this 
Fj,(t)|' acting can be expressed as 


{q-} + ['Sp® 1 { 4 ’} + {q’} 


Then it follows that 



(15) 


“ [‘P'lg {^'}g 


( 16 ) 
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The application of Eq. (6) now yields 

which removes the effect of the original branch from the inte- 
grated system. 


( 17 ) 


appendix b 

Rigid Arm Mass Branch 


Coupling 


Fig. IB 

Considering an offset mass 2 as shown in Fig. IB, the 


following definitions are assumed at 2 



Right hand rule is taken to define 6 , 0 and 6 
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[Mgl 


M, 


x2 


My2 


0 


M. 


z2 ' 

e-- 


( 


At 2, lumped mass In x, y and z direction are only taken, 
The kinetic energy of the mass 2 can be written as 

KE = I 


X, 


KE = ^ 


x2 

y2 

’.2 


But via the transformation 


M 


xl 


M. 


yi 


M. 


zl 


I 

r 


0 


0 


n . 

x„ 


'2 > 


x2 


y2 

z2 

J 


Cl) 


( 2 ) 


or 


A' 

X 


2 

^2 

z„ 


0 


> = 


x2 


y2 

z2 


0 

1 

0 

0 


0 I 0 
I 


0 

1 


-I 


I -I 

y X 


! 0 


! 0 


Cfig} = Etc] [6^] 


0 

1 

0 


3 

I 

i 

0 

0 

0 

1 


r> r\ 
^1 

^1 
z. 


0 


Xl 

®yi 

®zl 


(3) 


(4) 


REPRODUCII^i rn' OF 'O-F 
OKICI--: , 
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[tc] 






1 ^ 

0 

1 


1 


z 

y 

1 1 


0 


1 

Z 


X 

1 1 

i 


0 



X 


1 

1 

0 

^6"* 

0 ^ 
1 

0 

1 

0 

1 

0 

0 

1 


Substituting Eq. (4) in Eq. (1) yields 

= |([TC]{6^})*^fM2]([TC]t6^}) 

= I [6^}'^[MBC] {6^3 
where [MBC] - [tc]'^[M 2 ] [tc] 


(5) 


is the rigid arm mass branch coupling matrix, which describes the 
KE of mass 2 in terms of the velocity vectors of the degrees 
of freedom at point 1, 
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